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We present a measurement of D°-D° mixing parameters using the ratios of lifetimes extracted 
from samples of D° mesons decaying to K~tt + , K~ K + , and tt~-k + . Using 91 fb _1 of data collected 
by the BABAR detector at the PEP-II asymmetric-energy B Factory, we obtain a value Y = (0.8 ± 
0.4(stat.)l '|(syst.))%, which, in the limit of CP conservation, corresponds to the mixing parameter 
y = Ar/2r. Using the difference in lifetimes of D° and D° mesons, we obtain the CP-violation 
parameter AY = (-0.8 ± 0.6(stat.) ± 0.2(syst.))%. 

PACS numbers: 13.25.Ft, 12. 15. Ft, 11.30.Er 

To date there is no experimental evidence for mixing in the D°-D° system 0,0]- This is consistent with Stan- 
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dard Model expectations 0, Q , which correspond to a 
level of mixing beyond the reach of current experimen- 
tal precision. Among the more striking consequences of 
D°-D° mixing are different decay-time distributions for 
D° mesons that decay into final states of specific CP [f| . 
Measurable CP violation in D°-D° mixing would be ev- 
idence of physics beyond the Standard Model 0. 
The two D° mass eigenstates can be represented as 



|£>i> = p\D°) + q\D°) 
\D 2 ) = p\D°)-q\D°) , 



(1) 



where \p\ + \q\ = 1. It is traditional to quantify the size 
of D°-D° mixing in terms of the parameters x = Am/T 
and y = Ar/2r, where Am = mi — m 2 (Ar = Y\ — Y<i) 
is the difference in mass (width) of the states of Eq. 
and r = (Ti + ra)/2 is the average width. If cither x 
or y is non-zero, mixing will occur. The Standard Model 
expectation for the size of both is < 10 -3 

The effects of CP violation in D°-D° mixing can be 
parameterized in terms of the quantities 



and 



ff = arg 



qAj_ 

pA f 



(2) 



where A f = (f\H D \D°) (A f = (f\H D \D )) is the ampli- 
tude for D° (D°) decaying into a final state /. A value 
of r m ^ 1 would indicate CP violation in mixing. A 
non-zero value of iff would indicate CP violation in the 
interference of mixing and decay. Direct CP violation is 
expected to be small in the D°-D° system Q and is not 
considered here. 

D°-D mixing will alter the decay time distribution of 
D° and D° mesons that decay into final states of spe- 
cific CP. To a good approximation, these decay time 
distributions can be treated as exponential with effective 
lifetimes |7j 



[1 + r m (ycos(fi f - x sin tp f )]~ 



t = r° [l + r m 1 (y cos ipf +x sin </>/)] 1 , 



(3) 



where r° is the lifetime for the Cabibbo-favored decays 
D° -> K-TT+ and D° -> K+tt' and t+ (t") is the life- 
time for the Cabibbo-suppressed decays of the D° (D°) 
into CP-even final states (such as K~K + and ir~n + ). 
These effective lifetimes can be combined into the follow- 
ing quantities Y and AY: 



Y 



(r) 



1 



Ay 



A T 



(4) 



where (r) = (t + +t")/2 and A T = (t+ -t~)/(t + +t~). 
Both Y and AY are zero if there is no D°-D° mixing. 
Otherwise, in the limit of CP conservation in mixing, 
Y — y cos (p f and Ay = x sin ip / . 

We present a measurement of Y and AY obtained 
from a 91 fb _1 data sample collected on or near the 



T(45) resonance with the BABAR detector at the PEP-II 
asymmetric-energy e + e _ storage ring. 

The BABAR detector, a general-purpose, solenoidal, 
magnetic spectrometer, is described in more detail else- 
where 01 ■ Charged particles were detected and their mo- 
menta measured by a combination of a drift chamber 
(DCH) and silicon vertex tracker (SVT), both operating 
within a 1.5-T solenoidal magnetic field. A ring-imaging 
Cherenkov detector (DIRC) was used for charged-particle 
identification. 

Four independent samples of D° and D° mesons were 
used in this analysis. The first three samples (referred 
to as tagged) correspond to D mesons that decayed into 
K~ir + , K~K + , and ir~ir + |b| and include the decay 
D* + — » D°ir + to suppress backgrounds and distinguish 
D° from D°. These three samples were used to measure 
Y and AY. The fourth sample (referred to as untagged) 
consisted of K~K + decays that were not matched to a 
D* + decay and was used to measure Y. 

D° candidates were selected by searching for pairs of 
oppositely charged tracks of invariant mass near the ex- 
pected value for a D° meson. Each track was required to 
contain a minimum number of measurement points in the 
SVT and DCH. The two D°-candidate daughter tracks 
were fitted to a common vertex. The fit probability of 
this vertex fit was required to be larger than 1%. The 
interaction point (IP) was determined by calculating the 
point in space most consistent with the D° trajectory 
and the beam envelope (approximately 6 (im high and 
120 /im wide). 

Each D° daughter track was subjected to a likelihood- 
based particle identification algorithm. This algorithm 
relied on the measurement of the Cherenkov angle from 
the DIRC and on the energy loss (dE/dx) measured with 
the SVT and DCH. The identification efficiency was 
approximately 80% for tracks within the DIRC accep- 
tance with a tt^~ misidentification probability of about 
2%. The average tt^ identification efficiency was approx- 
imately 90%. 

To reduce combinatorial background that tended to 
accumulate at lower momenta, each D° candidate was 
required to have a momentum in the e + e~ center-of-mass 
frame greater than 2.4 GeV/c. This requirement was 
also effective at removing D° mesons originating from 
the decays of B mesons. 

The proper decay time and its measurement error a t 
for each D° candidate were calculated using the D° and 
IP vertex fits. The world average D° mass mo and 
the momentum of the D° were used to calculate the boost 
of the D° and to obtain the proper decay time. The dis- 
tribution of at, uncorrelated with true decay time, peaks 
at a value of 160 fs, and has a long upper tail. Poorly 
measured D° candidates with at > 500 fs (16% of each 
sample) were discarded. 

The decay D* + — > D°ir + is characterized by a tt + of 
low momentum (n s ). To increase acceptance, n s candi- 
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FIG. 1: The reconstructed D° mass distribution (points) su- 
perimposed on a projection of the mass fit (curve) for the 
four D° samples. The fit was performed within the restricted 
ranges of mass indicated by the vertical dotted lines. The 
portion of the sample assigned by the fit to the background 
is indicated by the shaded region. 



TABLE I: The four D° samples, their use, and, as calculated 
inside a ±15 MeV/c 2 mass window, their size, and purity. 



Sample 


Measures 


Size 


Purity (%) 


K~n+ 


r u 


265,152 


99.4 


K~K+ 


(t), A t 


26,084 


97.0 


7T~7r + 


(r), A T 


12,849 


87.9 


Untagged K~ K + 


(r) 


145,826 


68.1 



date tracks were not required to include DCH measure- 
ments. To improve momentum resolution, a vertex fit 
was used to constrain each n s candidate track to pass 
through the IP. If the fit probability of this vertex fit was 
less than 1%, the D* + candidate was discarded. 

The distribution of the difference in the reconstructed 
D* + and D° masses (5m) peaked near 145.4 MeV/c 2 . 
Backgrounds were suppressed by discarding D* + can- 
didates with a value of dm that deviated more than 1 
(2.5) MeV/c 2 from the peak for those ir s tracks measured 
with (without) the DCH. 

The D° mass distributions for the selected D° candi- 
dates are shown in Fig.^ Ample sidebands were included 
to measure the characteristics of the background. The 
peaks appearing above or below the D° mass were due 
to candidates with misidentified kaons or pions. For pre- 
sentation purposes only, we define those D° candidates 
with reconstructed masses within 15 MeV/c 2 of mr> as 
belonging to a mass signal window. The sizes and esti- 
mated purities of the four D° samples within this window 
arc listed in Table H 

An unbinned maximum-likelihood fit was used to ex- 
tract the lifetime from each D° sample. The likelihood 



function consisted of two decay-time distribution func- 
tions, one for signal and one for background. The sig- 
nal function was a convolution of an exponential and a 
resolution function that was the sum of three Gaussian 
distributions with zero mean. The widths of the first 
two Gaussians were proportional to a t whereas the width 
of the third, designed to describe mismeasurements, was 
not. The parameters in the fit associated with the signal 
for the K~ir + and untagged K~K + samples were the 
lifetime and the widths and relative proportions of the 
three Gaussians. The parameters for the tagged K~K + 
and 7r~7r + samples were the same except for the addition 
of A T . 

As in the signal likelihood function, the background 
function was a convolution of a resolution function and 
a lifetime distribution. The background lifetime distri- 
bution was the sum of an exponential distribution and a 
delta function at zero, the latter corresponding to prompt 
sources of background that originated at the IP. The 
resolution function consisted of the sum of four Gaus- 
sian distributions, the first three of which were similar to 
those of the signal. The fourth was given a fixed width of 
12 ps and accounted for a small number (< 10~ 3 ) of out- 
liers produced by long-lived particles or reconstruction 
errors. The additional fit parameters associated with the 
background included the fraction assigned to zero life- 
time sources, the background lifetime, and the relative 
size of the fourth Gaussian. 

To combine the signal and background likelihood func- 
tions, the reconstructed mass of each D° candidate was 
used to determine the probability that it was a signal D°. 
This calculation was based on a separate fit of the recon- 
structed D° mass distribution (Fig.Q. This fit included 
a resolution function composed of a Gaussian with an 
asymmetric tail designed to account for final-state pho- 
ton radiation. The mass fit for the tagged D° samples in- 
cluded a linear portion to describe the background. The 
slope of the background was constrained with D° can- 
didates in the 5m sideband (151 < 5m < 159 MeV/c 2 ). 
For the untagged K~ K + sample, the size of the radiative 
tail was taken from the tagged K~K + sample and the 
background was described by a quadratic function. 

The results of the lifetime fits are shown in Fig. [21 
Typical values for the fit parameters were a background 
lifetime similar to the D° lifetime and a third Gaus- 
sian width that was several times larger than the typical 
decay-time error. 

To ensure that the analysis was performed in an ob- 
jective manner, the values of the r , (r), and A T fit pa- 
rameters were hidden until the analysis method and sys- 
tematic uncertainties were finalized. 

Potential biases in Y and AY were investigated us- 
ing Monte Carlo (MC) samples produced by a GEANT4- 
based detector simulation and processed by the same 
reconstruction and analysis programs as the data. To es- 
timate the behavior of both signal and background, the 
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Decay Time (ps) 



FIG. 2: The decay-time distribution for the four D samples 
(points) within a ±15 MeV/c 2 mass signal window superim- 
posed on a projection of the lifetime fit (histogram). The 
shaded histogram is the portion of the sample assigned by 
the fit to the background. The points presented below the 
histograms are the difference between data and fit divided by 
the statistical error with error bars of unit length. 

equivalent of 70 fb _1 of continuum and T(45) MC data 
were studied. To augment these samples, separate MC 
samples were generated for specific decay modes of the 
D° meson. The values of Y and AY calculated from the 
MC samples (generated with Y = and Ay = 0) were 
consistent with zero within statistical errors. As can be 
seen in Fig.|I5J the fraction of prompt background is larger 
in the tt~tt + and untagged K~ K + samples than in the 
other two, a tendency that is accurately reproduced by 
the MC simulation. 

Potential inaccuracies in the simulation of tracking 
were explored by varying within current understanding 
the assumptions used in the MC simulation, including the 
center-of-mass boost and energy, the strength of the mag- 
netic field, and tracking resolution and efficiency. Small 
charge asymmetries (at the level of 2 to 4%) in the recon- 
struction efficiency of tt s tracks produced slightly differ- 
ent momentum and angular distributions for D° and D 
mesons. The influence of these effects on Y and AY" was 
checked by applying the lifetime fits on the data after 
weighting events to remove these asymmetries. 

The size and decay time characteristics of backgrounds 
in the data were determined in the likelihood fits with- 
out using any MC input. The MC samples were used, 
however, to determine how well the fits account for the 
presence of background. As part of these tests, the size 
of specific types of backgrounds was varied within uncer- 



TABLE II: Summary of systematic uncertainties. 





Change in Y (%) 


AY (%) 




Tagged 


Untag 


Tagged 


Category 


K~K+ 


TV 7T + 


K~K+ 


K~K + 7T"7T + 


Tracking 


±0.1 


±0.3 


±0.2 


±0.1 ±0.1 


Background 


+ 0.3 
- 0.5 


±0.5 


±0.3 


±0.2 ±0.2 


Alignment 


±0.1 


±0.1 


±0.1 


±0.1 ±°-° 


MC Statistics 


+ 0.4 
- 0.1 


+ 1.0 
- 0.1 


+ 0.4 
- 0.1 


< 0.1 < 0.1 


Quadrature Sum 


±0.5 


+ 1.2 
- 0.6 


+ 0.5 
- 0.4 


±0.2 ±0.2 



tainties in the MC sample by reweighting. In addition, 
MC parameters associated with the charm fragmentation 
function and final-state radiation in D° decay were var- 
ied. The resulting effect on the fitted lifetime is reported 
as one source of systematic uncertainty. 

Detector misalignment was another potential source of 
bias. Residual distortions of the SVT, even as small as a 
few microns, can produce significant variations in the ap- 
parent D° lifetime. Several studies were used to measure 
and characterize such distortions and strategies were de- 
veloped to correct them. One example was the study of 
proton tracks that were created by the interaction of off- 
energy beam particles and the beampipe. These tracks 
were used to measure the radius of the beampipe to a 
precision of a few microns, which limited the uncertainty 
in the radial scale of the SVT to 0.3%. 

Another example was a study of e + e~ — > e+e~ ± 
2(7r + 7r~) events in which the four pions were known to 
originate from the IP. By selecting oppositely charged 
pairs of these pions with opening angles similar to two- 
body D decays, it was possible to measure the apparent 
beam position as a function of £>°-candidate trajectory 
and calculate a correction to the D° lifetime. This type 
of correction nearly cancels in the lifetime ratio and in- 
troduces little systematic uncertainty in Y or AY. 

Because Y and Ay were measured from the ratio and 
asymmetry of lifetimes, systematic uncertainties from 
alignment that have a strong influence on r° did not make 
a large contribution to Y and Ay. Efforts to reduce these 
systematic uncertainties are still underway; therefore, a 
value of t° is not reported in this letter. A subsample of 
the data was used to verify that r° is consistent within 
uncertainties with the world average Q . 

The systematic uncertainties in Y and Ay are summa- 
rized in Tablelrfl The separate results for each sample are 
listed in Table ITTT1 with combined values that assume the 
same value of tpf for the K~ K + and 7r~7r + decay modes. 
All values are consistent with no mixing. Because it is 
derived from an asymmetry, the systematic uncertainty 
in Ay is considerably smaller than in Y . 

In summary, we have obtained a value of Y = (0.8 ± 
0.4(stat.)^Q^(syst.))% that is consistent with no mixing 
and is at least twice as precise as previous measurements 
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TABLE III: Summary of Y and AY results. The first error 
is statistical; the second, systematic. 



Sample 


Y (%) 


AY (%) 


K~K+ 


1.5 ±0.8 ±0.5 


-1.3 ±0.8 ±0.2 


7T~~7T + 


1.7 ± 1.2 tH 


0.3± 1.1±0.2 


Untagged K~ K + 


0.2 ±0.5 to.l 




Combined 


0.8 ±0.4 to.l 


-0.8 ±0.6 ±0.2 



of this type [H El El, all of which assumed CP 
conservation. We also obtain for the first time a mea- 
surement of Ay = (-0.8 ± 0.6(stat.) ± 0.2(syst.))%. 
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